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ABSTRACT 

We present an analysis of the neutral hydrogen (HI) properties of a fully cosmological 
hydrodynam ical dw arf galaxy, run with varying simulation parameters. As reported 
by iGovernat o et all (|2010!), the high resolution, high star formation density threshold 
version of this galaxy is the first simulation to result in the successful reproduction of 
a (dwarf) spiral galaxy without any associated stellar bulge. We have set out to com- 
pare in detail the HI distribution and kinematics of this simulated bulgeless disk with 
what is observed in a sample of nearby dwarfs. To do so, we extracted the radial gas 
density profiles, velocity dispersion (e.g. velocity ellipsoid, turbulence), and the power 
spectrum of structure within the cold interstellar medium from the simulations. The 
highest resolution dwarf, when using a high density star formation threshold compa- 
rable to densities of giant molecular clouds, possesses bulk characteristics consistent 
with those observed in nature, though the cold gas is not as radially extended as 
that observed in nearby dwarfs, resulting in somewhat excessive surface densities. The 
lines-of-sight velocity dispersion radial profiles have values that are in good agreement 
with observed dwarf galaxies, but due to the fact that only the streaming velocities of 
particles are tracked, a correction to include the thermal velocities can lead to profiles 
that are quite flat. The ISM power spectra of the simulations appear to possess more 
power on smaller spatial scales than that of the SMC. We conclude that unavoidable 
limitations remain due to the unresolved physics of star formation and feedback within 
pc-scale molecular clouds. 

Key words: galaxies: dwarf- galaxies: evolution - galaxies: formation - methods: 
N-body simulations 
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INTRODUCTION 



A tradit ional problem plaguing the simulation of disk galax- 
ies (e . g. iThacker fc Cou chman 200l! : ISommer-Larsen et al.l 



20031: 



Abadi et al 



2003; Govcrn ato et al.l 



2004 



20071: 



Robertson et al.l bOOJTlBailin et all |2005|; lOkamoto et all 



2005. : .Sanchez-Blazguez et al.ll2009l : IStinson et al.ll2010l . and 



references therein), within a cosmological context, has been 
the inability to recover successfully the properties of a truly 
"late-type" disk and, in particular, those with essentially no 
associated stellar bulge, similar to classical galaxies such as 
M33. 

Recent work by IGovernato et al.l l|20ld ). though, has 
produced what appears to be exactly such a bulgeless dwarf, 
via the imposition of a higher density threshold for star for- 



mation (100 cm~^, as opposed to 0.1 cm~^, as adopted in the 
aforementioned earlier generations of simulations) , and mass 
resolution that allows one to identify individual star form- 
ing regionslj The primary dwarf in their analysis (DGl) 
forms a shallow central dark matter profile and possesses a 
pure exponential stellar disk of radial scale r^^l kpc, with 
a stellar bulge-to-disk ratio _B/D~0.04 as determined from 
the i band light profile. 



^ The higher star formation density threshold can only be ap- 
plied because the high resolution of the simulation, coupled with 
heating from the UV background, ensures fragmentation does not 
occur at unresolved scales. 

^ In addition to the supplementary re-simulation of DGl de- 
scribed in § 2. 
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In what follows, we extend this work and examine in 
detail the cold neutral hydrogen (HI) gas content of the sim- 
ulated dwarf DGl along with its low star formation thresh- 
old analog, DGILT, and an updated version of DGl (called, 
nDGl) which employs high-temperature metal-line cooling 
and enhanced supernova energy feedback to compensate for 
the additional cooling. Our goal is simple: to determine if 
their HI gas properties agree with recent observational data 
to an equally successful degree as the stellar component. 
Studies such as the The HI Nearby Galaxy Survey (THINGS 
I Walter et al.l 12008;' ) provide excellent high resolution (spec- 
tral and spatial) data against which to compare our simula- 
tions. The gas properties of the simulations are compared 
directly wit h several of the most recent relevant empiri- 
cal datasets (IStaniinirovic et al.ll 19991 : iTamburro et al.ir2009l : 



lOIBrien^lalJboid), in order to ll^ ss both their stre ;;^s 
and weaknesses. 

The cold gas in galaxies is linked directly to under- 
lying star formation processes and associated interstellar 
medium (ISM) physics; any successful model of galaxy for- 
mation should adopt a holistic approach, examining both 
the gcis and star properties in consort. We describe the ba- 
sic properties of our simulations, before detailing the anal- 
yses undertaken; we will present results pertaining to the 
radial distribution of cold gas within the disks associated 
with DGl, DGILT, and nDGl, spatially-resolved velocity 
dispersion maps of the cold gas, and the spatial distribu- 
tion of power encoded within the structure of the ISM. We 
end with a summary of our findings, discussing both the 
strengths and weaknesses of the current simulations. 



2 METHOD 



2.1 Simulations 



We have made use of the recent iGovernato et al.l (|201G| ) 
simulations which produced, for the first time, a late-type 
dwarf spiral with no associated stellar bulge. A full de- 
scription of the simula tions' characteristics is provided by 
IGovernato et al.l l|2010r ). but for context, it is useful to sum- 
marise their primary traits. 

Using the N-body-|-SPH (|Monaghanlll992l ') code GASO- 
LINE, a low resolution (25 Mpc box, sufficient to provide 
realistic torques for these dwarfs), dark matter only sim- 
ulation was used to identify 3.5 xlO'^" M© (virial) halos 
(with typical spin paramters A=0.05) for potential (high 
resolution) re-simulation using a volume renormalisation 
technique (i.e., "zoom" simulation). New initial conditions 
were then re-constructed for the primary target halo (called 
"DGl"), using the relevant low-frequency waves associated 
with tidal torquing in the low resolution "parent" simula- 
tion, but now enhanced with higher spatial frequencies gen- 
erated after tracing the present-day particles back to the 
relevant Lagrangian sub-region within the parent. The mass 
distribution was then sampled at higher resolution in the 
regions of interest, and more coarsely, further away from 
the identified halo. Both DGl and nDGl have a force res- 
olution of 86pc, while that of DGILT is somewhat lower 
(116pc); the initial baryonic (dark) particle mass for DGl 
and nDGl is 3300 M© (16000 Mq), while for DGILT it is 
7800 M0 (37000 M©). At z=0, the i-band luminosities of 



DGl, nDGl, and DGILT are M,=-16.5, -15.8, and -19.1, 
respectively. 

We should re-iterate that each of the three simulations 
described here (DGl, DGILT, and nDGl) use the same dark 
matter halo / assembly history, and differ primarily only in 
their treatment of the baryonic physics associated with star 
formation - i.e., either supernova energy feedback efficiency 
(DGl vs nDGl) or star formation density threshold reso- 
lution (DGl vs DGILT). DGl was simulated using a star 
formation density threshold of 100 cm~^, typical of the den- 
sities encountered in giant molecular clouds, rather than the 
canonical value adopted in earlier simulations (0.1 cm~'')|j 
Other than the increased density threshold, two additional 
parameterisations were adopted, within the context of the 
feedback formalism employed: the star formation efficiency 
(eSF=0.1) and the fraction of supernova (SN) energy cou- 
pled to the ISM (eSN=0.4). Th e star formation and f eedback 
are modelled as described in IStinson et al.l (|2009l ). With- 
out any additional ad hoc adjustments, this high density 
threshold led to bulgeless dwarf spirals (akin to the clas- 
sic prototype, M33) with flat (non-centrally concentrated) 
rotation curves (again, for the first time). Alongside our 
analysis of the high-threshold DGl simulation, we provide 
a parallel analysis of two other simulated dwarfs, DGILT 
(the lower-threshold analog, which uses the aforementioned 
canonical 0.1 cm"'' threshold, and a star formation effiency 
eSF=0.05, with the same initial conditions as that used for 
DGl), and an updated version of DGl, nDGl (again with 
the same initial conditions as DGl and high density thresh- 
old of 100 cm ~^, but now with high-temperature metal-line 
cooling, after IShen et al.l (|2010|), and increased thermal en- 
ergy coupling to the ISM (eSN=l)), in order to better assess 
the role played by star formation threshold and feedback in 
"setting" the gas properties of the respective simulations. 

To foreshadow the discussion which will follow, perhaps 
the most problematic aspect of the current analysis is the 
uncertain numerical "leap-of-faith" that must be made in 
associating the typically 7000— 8000K SPH gas particles, re- 
gardless of their local density (~0. 1 — 100 cm~^), with star 
formation (which in nature occurs in clouds and cores with 
temperatures 2—3 orders-of-magnitude lower than this). Un- 
til the effects of cooling by molecular hydrogen are incorpo- 
rated fully within Gasoline, this remains an unavoidable 
limitation of our modeling, but fortunately one whose effects 
are known and well-understood. We return to this point in 
§ 2.2 and § 3.3. 



2.2 Analysis 

The cold gas properties of DGl, DGILT, and nDGl are 
compared directly with those from comp arable dwarfs in 
The HI Nearby Galaxy Survey (THI NGS: ITamburro etlo] 
|2009|). in addition to the sam ples of lO'Brien et al.1 (l2010h 
and IStanimirovic et al.l l|l999r ). The bulk properties of DGl 
(e.g., mass, luminosity, and gas fract ion) are con s istent 
with those observed in nature (e.g. IWalter et al.l l2008l : 



^ Gasoline employs an ideal gas law equation of state 
llWadsIev et al.l l2004). and the mea n molecular weig ht is implic- 
itly solved for and allowed to vary llShen et al.ll201(]|) . 
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Ivan den Bosch et alj 1200 ih . and its present-day star for- 
mation rate (~0.005 Mo/yr) and luminosity (Mi~— 16) 
are (specifically and directly) comparab le to those of the 
three dwarfs from lTamburro et al.1 ((20091), with Holmberg II 
(HoII) being perhaps the closest direct analog (and, as such, 
being the empirical counterpart to which we will refer DGl 
most often). As noted earlier, the properties which we de- 
rive include the radial extent, the velocity dispersion as a 
function of galacto-centric radius, and the power spectrum 
of the ISM. 

In our work, unless otherwise stated, we label "cold 
gas" those SPH pa rticles with temper atures less than 
Tmax=15000 K f after IStinson et al.1 (|2006l )). The bulk of the 
gas in DGl (uDGl) lies near 7000K (9000K), which at face 
value would appear to be more appropriate for the warm HI 
phase of the ISM, rather than the cold, star-forming, gas, to 
which we have associated star formation within the simula- 
tion. However, our cooling, despite the inclusion of metal- 
line cooling, is limited primarily to hydrogen and helium 
cooling, which can only cool gas down to these temperatures, 
and as emphasised in Stinson et al (2006; §5.1.1), we are 
averaging over scales much larger than individual star form- 
ing cores. The effect of varying this maximum temperature 
th reshold (Tmax) for st ar formation was examined in detail 
bv lStinson et al.l (|2006l ). to which the reader is referred. We 
can summarise that analysis by stating that provided Tmax 
is chosen to be not too similar in value to that of the mean 
temperature of the gas particles, its sp ecific value does not 
critically affect star formation (see also. IShen et al.l (2010 )'). 
Efforts are underway to implement molecular hydrogen cool- 
ing within Gasoline, after which a quantitative comparison 
with our results can be undertaken. 

DGILT, the low density threshold analog to DGl, is 
analysed in parallel, to provide s omething of a canonica l 
"control" sample. As described in iGovernato et al.l (|201Cl ). 
the properties of DGILT (e.g., rotation curve, dark mat- 
ter density profile, bulge-to-disc ratio) are not well-matched 
to those observed in nature, due to the traditional lim- 
itations that the new suite of simulations were designed 
to overcome in the first place. As a juxtaposition to DGl 
though, it is invaluable. The present day star formation rate 
(0.2 Mo/yr) and luminosity (Mi=— 19.1) are much higher 
than that of DGl (and the associated stellar mass is corre- 
spondingly a fac t or of ten higher), driven (as described by 
IGovernato et al.l (|2010|)) by its adoption of the lower star 
formation threshold (see Fig 1). 

For our analysis, we have generated a new variant of 
DGl (labelled nDGl), employing both the same initial con- 
ditions and the higher star formation threshold (100 cm~^). 
As alluded to earlier, where nDGl differs from its prede- 
cessors is in its inclusion of metal-line cooling (following 
IShen et al.l (12010|)) and a more efficient coupling of SN ther- 
mal energy to the ISM; qualitatively, we can anticipate this 
leading to a somewhat more turbulent ISM. On the whole, 
the star formation rate of nDGl is suppressed relative to 
DGl, but extends to lower redshifts (Fig 1, where one can 
see that the star formation rate from 8<i<10 Gyr is ~10x 
higher in nDGl than in DGl); its luminosity is, not sur- 
prisingly, somewhat lower than that of DGl (Mi = — 15.8, as 
opposed to Mi=— 16.5), considering its stellar mass is a fac- 
tor of two lower (M,«2.1xl0* Mq vs M.Ri4.4xl0* Mq). 

Zeroth (density), first (velocity), and second (velocity 




Time [Gy 



Figure 1. The star formation rates of nDGl (solid line), DGl 
(dot-dashed line), and DGILT (dashed line). Star formation in 
nDGl is suppressed overall, relative to DGl, but extends ~2 Gyrs 
beyond the cessation of bulk star formation in DGl (in the 
range S^t^lO Gyrs). There is intermittent star formation in both 
dwarfs up to the present day. but it has been consistently low for 
the past ^S Gyrs in nDGl and ~5 Gyrs in DGl. The star forma- 
tion history of DGILT is overall considerably higher than its two 
higher density threshold analogs. 



dispersion) moment maps of the simulated neutral hydrogen 
distributions were generated using tipsMj , after matching 
the ~ 40° inclination of the dwarfs from the lTamburro et al.l 
(|2009l ) THINGS sample (which, again, includes HoII, our 
primary analog against which our simulations will be com- 
pared, as noted in § 2). The conversion from "cold gas" 
to "HI" within GASOLINE suffices for the purposes out- 
lined here; the values derived are close to the values one 
would predict under the assumption of combined photo- and 
coUisional-ionisation equilibrium. All our results were cross- 
checked using both cold gas and HI moment maps, in ad- 
dition to further cross-checks undertaken after eliminating 
high column density HI gas for which the conversion from 
cold gas to HI is most insecure. The results described here 
are robust to these choices, and for expediancy are not dis- 
cussed further. 

Our velocity dispersion analysis made use of the second 
HI moment map derived from the line-of-sight dispersion 
map produced from viewing the DGl, nDGl and DGILT 
simulations with an inclination angle matching that of HoII. 
For the analysis of the distribution of structural "power" 
within the cold ISM of the simulations, we again used the 
zeroth HI moment maps and their Fourier Transforms, and 
compared the i nferred power law spectra with that derived 
for the SMC bv lStanimirovic et al.1 (|l999l ). 



3 RESULTS 

3.1 Radial Density Profiles 

We first confirmed independently that the stellar light 
associated with DGl was indeed consistent with a pure 

www-hpcc . astro .Washington. edu/tools/tipsy/tipsy .html 
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exponential of scalelength '^l kpc (i.e. bulgeless) disk 
(akin to the Typ e I pr ofiles categorised by, for example, 
IPohlen fc Truiilld (|200a)); as shown in the lower panel of 
Fig 2, this was the case. DGfLT also has a radial (stellar) 
scalelength of ~1 kpc, but shows the classical "problem" 
of possessing an substantive stellar bulge within the inner 
kpc (_B/Z)~0.2). The stellar disk component of nDGl is 
not well-represented by a single pure exponential (cf. DGl); 
instead, its surface density profile shows a deficit of mat- 
ter (and light) in the outskirts of the stellar disk (beyond 
a so-called "break radius" at ~2— 3 kpc), consistent with 
the more common Typ e II profiles observed in nature (e.g. 
iPohlen fc Truiilld l2006l : ISanchez-Blazguez et al]|2009l ): the 
inner and outer parts of the nDGl stellar disk show ra- 
dial scalelengths of ~2 kpc and '^l kpc, respectively. The 
bulge-to-disc ratio of nDGl matches formally that of DGl, 
although it is also readily apparent that the surface den- 
sity (and light) profile of nDGl shows a high-density stel- 
lar "core", in which -lO'^ M© (~10% of the nDGl steUar 
mass, as a whole) is concentrated within the inner 100 pc. 
Importantly, this stellar "core" is inconsistent with a bulge. 
Instead, it consists of a large cluster of stars that was formed 
in the disk during a merger at high-redshift, and traveled in- 
ward with time so that at z=0 it is close to, but not located 
at, the dynamical center of the galaxy (i.e., it can be seen 
to rotate about the galaxy center). 

The cold gas of DGl displays a rapid increase in den- 
sity within ~lkpc. Exterior to this is an extended disk with 
an exponential scalelength ra ~6 kpc; the cold gas disk 
truncates at ~lrd, so mew hat short of those o bserved by 
iTamburro et al.l (|2009l 'l and lO'Brien et all (|2010(1. where the 
respect ive HI disks are traced out to ~2— 6 r^. iBigiel et al.l 
Il2008a|) showed that there is an empirical HI upper limit 
encountered in nature - J^hi ^9 Mo/pc^. This upper limit 
is represented by the horizontal line in the upper panel of 
Figure 2. Because we do not yet resolve the microphysics 
associated with molecular processes on parsec-scales, one 
might ascribe some fraction of the cold gas in the simu- 
lation (part i cularly that above the upper limit observed by 
iBigiel et al.l (l2008a|).) to molecular gas. Again using the re- 
sults from Bigiel et al. for the fraction of H2/HI as a function 
of radius (see their Figure 13), we can verify that the high 
density gas interior to 1 kpc is consistent with being molec- 
ular gas. In particular, r25, the isophotal radius correspond- 
ing to 25 mag/arcsec'^ , is 2.0 kpc for DGl. Assuming that 
as much of the gas can be ascribed to HI as possible (i.e., 
the upper limit of ^hi = 9 Mq/pc^), then the results from 
Bigiel et al. suggest that 7.9 Mq/pc^ would typically be in 
molecular gas at the innermost radius of DGl, dropping to 
7.2 Mo/pc^ at 0.8 kpc, and declining radidly to <0.1 H2/HI 
at 2 kpc. That is, while the total amount of gas in DGl is 
consistent with empirical bulk scaling relations, and the gas 
within 1 kpc is consistent with being mostly molecular, the 
cold gas surface densities beyond 1 kpc are too high relative 
to naturelj 

It is difficult to interpret the source of this excess gas. 



^ We note that while this is a sample of one, additional sim- 
ulatio ns from the same suite (e.g., DG2 from I IGovernato et al.l 
[201C|)) show the same behaviour; wc have chosen to focus only 
upon DGl, for clarity. 



Perhaps this is gas that should instead be lost from the 
galaxy in winds? While it may be tempting to suggest that 
this gas is overly concentrated, comparison of the cold gas 
scale lengths for these simulated galaxies (which has been 
fit beyond r25) to the scalelengths beyond r25, in the sample 
of I Bigiel et al.l l|2010l ). suggests that the excess gas in these 
simulations is actually too extended compared to real galax- 
ies. Alternatively, as discussed below for the case of DGILT, 
additional star formation in the outskirts of the simulated 
galaxy disks could decrease th e surface densi t y of g as (as it 
goes instead into stars). While iBrooks et al.l (201l|) showed 
that the _B-band scale length of DGl is comparable to ob- 
served dwarf galaxies, a factor of 1.5 to 2 increase in size 
is still allowable to be fully consistent with nature. In fact, 
preliminary tests of molecular cooling and star formation in 
Gasoline suggest that the star formation is more extended 
at 2=0. Hence, the addition of H2 to these simulations may 
alleviate the problem of this excess gas. 

As was the case for the stellar light, the disk of nDGl is 
better represented by a "broken", or two-component expo- 
nential, with inner and outer disk scalelengths of ~2 kpc and 
~1 kpc, respectively (with the break occurring near a galac- 
tocentric radius of ~3 kpc). The arguments of the previous 
paragraph concerning the excess surface density of cold gas 
in DGl applies obviously to nDGl, as well. 

Conversely, the cold gas in the disk of DGILT extends 
radially to ~8 kpc with an essentially flat density profile 
(formally, with a radial scalelength of ~18 kpc - i.e. , the 
gas disk truncates near ~0.5rd - again, short of the typical 
disc in nature, but since the profile is so fiat, the formal ex- 
ponential "scalelength" is somewhat ill-defined). Like DGl, 
DGILT also shows a high density cold gas "core" (of mass 
~2xl0® M0), although it is somewhat more extreme, in 
the sense of it being concentrated solely within the inner 
~100 pc (note that this is within twice the force softening 
length) . Being more extended, and the gas fraction being an 
order-of-magnitude lower (Gwernato ct al. (2010); Tbl 2), 
it is not surprising that the cold gas surface density profile 
of DGILT is consistently a fac tor of ~3 x lower th an the em- 
pirical upper limit derived by (JBigiel et al.lbOOSg 'l. However, 
this result should not be interpreted to mean that DGILT 
is the more real i stic v ersion of this g a laxy simulation. As 
iGovernato et all l|2010t ) and lOh et al.l (|201ll ') have demon- 
strated clearly, the mass of this galaxy is overly concen- 
trated, with a large bulge and peaked inner rotation curve 
that are inconsistent with observed galaxies in the same 
mass range. 



3.2 Velocity Dispersion 

We next undertook an examination of the velocity dis- 
persion of the HI disks of DGl, nDGl and DGILT, to 
make a comparison with those o bserved in various sam- 
ples of dwarfs in the literature (ICrosthwaite et al.l |2000| . 
I2OOII : ITamburro et all l2009l : lO'Brien et al.ll2010l ). Observa- 
tions show that independent of present-day star formation 
rate, luminosity, or mass, disks possess a characteristic veloc- 
ity dispersion of ~8— 10 km/s, rising to ~12— 15 km/s in the 
inner star-forming regions (i.e. within r25, the isophotal ra- 
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O DG1 Cold Gas (r^ = 5.7kpc) 

XnDGI Cold Cos (r,, = 2-3kpc) 

X nDGl Cold Gas (r„„ = 0.9kpc) 

ADG1LT Cold Gas {r^' = lB.3kpc) 




Figure 2. Radial gas (top) and stellar (bottom) density profiles 
for the simulated dwarfs DGl (diamonds), DGILT (triangles), 
and nDGl (crosses). The thick overplotted lines show the expo- 
nential fits to the distributions, from which the noted scalelengths 
were derived. The stellar component of DGl obeys a pure expo- 
nential of scalelength ~1 kpc, with no evidence for a central bulge, 
while both nDGl and DGILT show central cores. Both the stellar 
and cold gas components of nDGl are best represented by double 
exponentials, with a break between the two near ~3 kpc. The 
cold gas of DGl is distributed in a more extended exponential 
disk component of scalelength ~6 kpc, while that of DGILT is 
~18 kpc. The horizontal line in the upper panel corresponds to 
the empirical u pper limit to HI en countered in nature, from the 
THINGS work l|Bigiel et al.ll2008a^ . 



dius corresponding to 25 niag/arcsec^)|j A typical radial ve- 
locity dispersion distribution is sho wn in Fig 3 for HoII (p lus 
signs), from the THINGS sample (|Tamburro et al.ll2009h . 

In addition to the curve for HoII, in Figure 3 we also 
show the corresponding velocity dispersion profiles (line-of- 



^ At the resolutions at which we are working (~100 pc), the ve- 
locity dispersions of the molecular and neutral gas are not dra- 
matically different - /citepCrosOO,Cros01. 



sight, assuming again a ~40° inclination, similar to that 
of HoII) for DGl (open diamonds), nDGl (crosses), and 
DGILT (triangles), derived from the SPH gas particles' 
s tream ing velocities (see below, and Ivan den Bosch et al.l 
(2002)), and for DGl (filled diamonds), taking into account 
said particles' thermal velocities. Circular annullj projected 
on the inclined galaxy were used to set the bins. 

For typical Milky Way-scale simulations, the thermal 
broadening component is often neglected, since the 'stream- 
ing velocity' of the SPH particle usually dominates over 
the 'thermal component'. For our simulated dwarfs, this is 
clearly inadequate, as the streaming velocity dispersion can 
be much smaller than the relevant thermal velocity disper- 
sion. To incorporate the latter, we f ollow the procedure out- 
lined bv lvan den Bosch et al.l l|2002r ) (§2.3) and note that the 
velocity of each particle can be written as v = u + w, where 
u is the mean streaming velocity at the location x and w is 
the particle's random (thermal) velocity. Because SPH only 
tracks the streaming motions of the particles, we make use 
of the internal energy of each particle, in order to derive an 
appropriate random component to apply to each particle. In 
practice, we draw random velocities for each Cartes ian coor- 
dinate from a Gaussian of dispersion a = ^kT j \i, and add 
those to each of the coordinates of the streaming motion, 
where T is the temperature of the gas particle (typically, 
~7000— 9000 K, for our simulations), k is Boltzmann's con- 
stant, and /i is the mean molecular weight of the gas. 

Without the inclusion of thermal broadening, both DGl 
and nDGl show extremely (and unphysically) kinematically 
cold interstellar media compared to DGILT and, more im- 
portantly, dwarfs in nature (compare the crosses and open 
squares of Figure 3 (simulations) with those of the plus sym- 
bols (observations) for a graphic example of the mismatch 
between unphysical streaming velocity dispersions and those 
encountered in nature). This is not to imply, however, that 
DGILT as presented in Figure 3 is physical. First, and most 
importantly, as a heady not ed in § 3. 1 and, especially, by 
ICovernato et al.1 (|20ld ) and lOh et ai] (l2011r ). the rotation 
curve and dynamics of DGILT are problematic, as is the 
associated significant overproduction of the stellar bulge. 
As can be seen in Fig 1, DGILT has a star formation rate 
two orders of magnitude larger than DGl or nDGl; while 
this does not impact upon its consistency with the stellar 
mass-metallicity, luminosity-metallicity, or HI gas fraction- 
luminosity scaling relations, it does worsen significantly the 
consistency with the dynamica l-to-stellar mass ratio distri- 
bution of iBlanton et al.l (|2008l ). This large star formation 
rate drives more turbulence, leading to the large streaming 
velocities for this simulation. We have not included the ther- 
mal component for DGILT in Fig 3, as doing so would only 
increase its velocity dispersion from ~12 km/s to ~14 km/s. 
The inferred line-of-sight velocity dispersion profile for DGl, 
after application of the above thermal broadening (which ef- 
fectively amounts to a (t~7— 9 km/s broadening of the essen- 



^ Technically, elliptical annuli should be used, but our results are 
not sensitive to this choice, at these inclination angles; in addition, 
we re-measured the velocity dispersion profile on the raw THINGS 
data for HoII using circular annuli, to ensure self-consistency with 
our analysis of the simulations. 
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tially negligible ~1 km/s streaming motions), is represented 
by the filled squares in Figure 3. 

The characteristic velocity dispersions of the cold gas 
within DGf and nDGl a re comparable to those encountered 
in nature (~8— f km/s - iTamburro et alj (|2009')) when ther- 
mal velocities are considered. The thermally broadened ve- 
locity dispersion profile of DGf shows a few enhanced fea- 
tures (near 0.5r25). These are due to high temperature gas 
particles in and around superbubbles blown by SNe feedback 
(discussed further below and shown in Figure 4). By de- 
sign, including a random thermal component to the velocity 
dispersion accentuates these features. However, by chance, 
the particular timestep we examine here for uDGf does not 
show any bubbles (though does at previous timesteps), and 
hence no thermal features are introduced into the profile of 
this simulation. As can be seen from the streaming-only pro- 
files for these galaxies, both have slightly higher macroscopic 
velocity dispersions in the inner few hundred parsecs. How- 
ever, in DGf this gas is ~35% hotter than the rest of the 
disk, while in uDGf it is cooler by a similar factor. Figure 3 
shows that, when this is considered in the thermally broad- 
ened velocity dispersions, it has the effect of maintaining 
the higher velocity dispersion structure in the inner region 
of DGf, while "washing out" the inner structure in uDGf . 
This result highlights a conundrum in terms of comparing 
the velocity dispersion profiles of these dwarf galaxy simu- 
lations to real dwarfs. 

A more subtle effect of imposing the random thermal 
velocity perturbation to each particle's streaming motion is 
that the velocity ellipsoid of the cold gas becomes neces- 
sarily isotropic, disguising any anisotropies that might have 
been present in the streaming motions (i.e., young stars, and 
the cold gas from which they formed, will necessarily have 
different velocity ellipsoids). For example, for DGf (nDGl), 
the radial, azimuthal, and vertical velocity dispersions in- 
ferred from the cold gas particles' streaming motions, mea- 
sured at ^0.5rd, are (Tr-~4 km/s (^^6 km/s), CT0~3 km/s 
(~6 km/s), and a^^l km/s (~2 km/s) -i.e., (Tr:cr^:o"z«3:3:f 
(anisotropic). After thermal broadening, the derived re- 
spective velocity dispersions are CTr~8.5 km/s (~fO km/s), 
cr0«8 km/s (~fO km/s), and CTz~7.5 km/s (~8.5 km/s) - 
i.e., (Tr-:o'</>:crz~f :l:f (isotropic). What this means is that an 
unavoidable outcome of our current inability to resolve pc- 
scale molecular heating and cooling processes within the 
simulations is the lack of any significant correlation between 
velocity dispersion and galactocentric radius and/or under- 
lying star formation. Until we can resolve densities (and tem- 
peratures) corresponding to the cores of molecular clouds, 
this apparent mismatch between observations and simula- 
tions would appear difficult to avoidlj 



* It might be tempting to conclude that since the enhanced feed- 
back did not result in a significantly higher linc-of-sig ht velocity 
disper sion, this is consistent with the earlier work of iDib et al.l 
(I2OO6I ) and iPetric k, RupenI |2003), who concluded that super- 
nova feedback alone was insufficient to provide turbulent heating 
to the cold ISM in excess of a few km/s; in light of the fact that 
we are not resolving the ISM heating and cooling processes at pc 
and sub-pc scales, we feel it premature to draw such a conclusion 
from this aspect of our analysis. 



n DG1 w/out broodening 
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Figure 3. Radial behaviour (in units of the B-band r25 - i.e. the 
isophotal radius corresponding to 25 mag/arcsec^ or, roughly, to 
the extent of the star forming disk) of the HI line-of-sight velocity 
dispersion of the DGl (open squares), DGILT (open triangles), 
and nDGl (crosses) simulations, d erived from the SPH gas pa r- 
ticles' 'streaming velocities' (after Ivan den Bosch et al.l (12003)), 
in addition to the true HI line-of-sight velocity dispersion profile 
for DGl (filled squares) and nDGl (plus signs), after correcting 
the streaming velocities isotropically for their internal thermal 
energies. Also shown is a representative dwarf spiral from the 
THINGS (Tamburro et al. 2009) sample (HoII: open diamonds), 
note: r25 is 2.0 kpc, 5.5 kpc, 1.4 kpc, and 3.3 kpc, respectively, 
for DGl, DGILT, nDGl, and HoII. 



3.3 Povirer Spectrum and Superbubbles 

Following IStanimirovic et al.l l|f993 ). we generated the 
Fourier Transform of the HI moment zero maps of DGl, 
uDGf, and DGILT - each shown in Figure [4] at the same 
spatial scale (14x14 kpc) with the same limiting HI col- 
umn density (N(HI)>f xfO^'' cm~^) - after first convolv- 
ing the maps with a 100 pc Gaussian, to mimic the typi- 
cal beam-smearing pre sent within THINGS data for HoII 
(|Tamburro et alJ I2OO9I ) . Circular annuli in Fourier space 
were then employed to derive the average power in the struc- 
ture of the ISM on different spatial scales. Figure[5]shows the 
derived power spectra for the simulations DGl, nDGl, and 
DGILT, and that for the Small Magellanic Cloud (SMC), 
re-derived for self- consistency, using the HI datacube of 
IStanimirovic et al.l (|1999. ). Grossly speaking, the distribu- 
tions can be represented by a power law of the form Pock'' , 
with 7=-3.5 for DGl, 7=-3.4 for DGILT, and 7=-4.2 for 
nDGl, and 7 =— 3.2 for the SMC (cons istent with that found 
originally bv IStanimirovic et al.l (Il999l ). and consistent with 
the power spectrum expected when HI density fluctuations 
dominate the ISM structure, rather than turbulent velocity 
fluctuations, which dominate the spectrum when isolating 
'thin' velocity slices). 

There are several points to highlight from Fig [5] (i) the 
SMC shows no evidence for departure from a pure power 
law, and hence there does not appear to be any obvious pre- 
ferred HI cloud size in nature; (ii) broadly speaking, both 
DGl and DGILT are shallower than nDGl (i.e., possess 
more power on smaller scales, rather than larger, relatively 
speaking); put another way, the enhanced feedback associ- 
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Figure 4. Neutral hydrogen (HI) rtioinent zero maps of the three simulations analysed here - from left to right; DGILT, DGl, and 
nDGl. Each panel has dimensions 14x14 kpc; a lower column density threshold of N(HI)=lxlO^^ cm~^ was employed for each map. 



ated with nDGl shifts power in the simulated ISM from 
smaller scales to larger scales, just as one might expect; (iii) 
each of the simulations shows a greater departure from a 
pure power law, than does the SMC; the most obvious depar- 
ture from a power law is perhaps seen in the enhanced power 
on ^--^400— 500 pc scales in nDGl. This enhanced power cor- 
responds to the "radial cadence" , or frequency, of the tighly- 
wound spiral structure in the inner few kpcs of the simula- 
tion (apparent in the right-most panel of Figure 4). 

Finally, from the present-day moment zero column den- 
sity map of DGl (middle panel of Figure 4), we identified 
13 SNe-driven superbubbles in its cold ISM. While we do 
not wish to belabour the point when employing such small- 
number statistics, it is re-assuring to note that upon plotting 
the superbubble size distribution, the data were consistent 
with a power law slope between —1.5 and —2.0 (dependent 
upon normalisation). Such slopes are entirely consisten t with 
those observed in nearby dwarfs (|Oev fc Clarkelll997l ). 



4 DISCUSSION 

One immediate concern arising from our analysis relates to 
the issue of extracting "neutral hydrogen" from the simu- 
lations' "cold gas" (which in some sense consists of both 
molecular and neutral hydrogen). Because the high-density 
regions within the simulation have densities more akin to 
molecular, rather than neutral, clouds, it is important to 
explore the definition of "neutral" employed here|f| To do 
this, we re-generated HI moment maps, but now restrict- 
ing the gas included to only those particles with densities 
near the classical value of ~0.1 cm~^. As expected, this 
eliminated the unrealistically high neutral hydrogen col- 
umn densities in the highest density regions, but at the 
expense of leading to vertical density profiles that bore lit- 
tle resemblance to th e Gaussian profiles observed in nature 
(|0'Brien et al.ll20ld ). Such an extreme "cut" to the defi- 
nition of neutral hydrogen also led to a radial profile that 



^ In large part, this was motivated by the fact that in "column 
density space" , these high-density regions possess column densi- 
ties close to 10^^ cm~'^, higher than those observed in nature; 
this is a limitation of the conversion employed within GASOLINE. 



A DGILT 7 = -3.4 
ODGl 7=-3.5 
7=-3,2 
7 = -4-2 




450 
Spotiol Scale (pc) 



Figure 5. Spatial power spectra of the cold ISM of DGl (dia- 
monds), DGILT (triangles), nDGl (crosses), and the SMC (plus 
signs). Power law slopes of —3.5, —3.4, —4.2 and —3.2 are over- 
plotted for DGl, DGILT, nDGl, and the SMC, respectively. The 
"brea k" in the SMC powe r spectrum is due to a missing baseline 
in the lStanimirovic et al.l ||l999) ATCA dataset. The power spec- 
tra for the three simulations have been truncated at ~2 resolution 
elements (2*FWHM of the adopted Gaussian beam: ~200 pc). 



bore little resemblance to an exponential. We found no den- 
sity cut which impacted favourably on the observable prop- 
erties of DGl. For these simulations, because density and 
temperature are closely correlated in the relevant regime 
(T<30000 K; p>0.001 cm"^), the above analysis is degen- 
erate to cuts in volume density or temperature. 

It is important to note that the primary process re- 
sponsible for driving bulk properties in the simulatio n is the 
star f ormation and feedback prescription. Governat o et al.l 
(|20ld ) demonstrated that star formation had a larger effect 
on the rotation curve of our simulated galaxy than resolu- 
tion (see their Figure 5). The gas properties presented in 
this paper are primarily the result of the star formation pre- 
scription, and thus it is imperative to use a star formation 
and feedback prescription that is physically motivated. Un- 
til metal-dependent H2 creation and cooling is added to the 
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simulations, it is not clear how much HI, as opposed to H2, 
should be present in the simulation, how it might be dis- 
tributed as a function of radius, and what impact it will 
have on the resulting disk. 

After applying the physically-motivated ~8 km/s ther- 
mal broadening to the Cartesian coordinates of the SPH 
particles' streaming motion, the inferred characteristic ve- 
locity dispersions for the cold gas were a reasonable match to 
those observed in nature (albeit, at the unavoidable expense 
of recovering any correlation between velocity dispersion and 
galactocentric radius and/or global star formation in the 
disk, in addition to the imposition of an isotropic velocity el- 
lipsoid to the cold gas, and the young stars which form from 
this gas) . Beyond the aforementioned issue of the lack of a 
self-consistent treatment of molecular cooling processings on 
sub-parsec scales, one must also be aware that at the reso- 
lutions of these simulations, we are still missing unresolved 
star forming regions and associated turbulence. The nature 
of these missing sources is an active area of debate, but 
magnetorotational instability (MRI) is one of the favoured 
mechanisms capab le of providing a non-negligible amount 
of turbulence fe.g. IWang fc AbelllJOOgl : iPiontek fc Ostrikeij 
l2007l : lMac Lowlboog ) 

Enhancing the supernovae energy feedback, as was done 
for simulation nDGl, at these resolutions, had a marginal 
impact on the SPH particles' streaming velocities (at the 
~20% level), which in turn meant that its impact on the ve- 
locity dispersion profiles was also minimal. This is not sur- 
prising, as the increased energy deposition was used in order 
to offset the effect of the newly included high-temperature 
metal-line cooling. Without the inclusion of extra SN energy, 
the additional cooling that comes from metal lines leads to 
more star forma tion than in the case of DGl. As shown by 
lOh et al.l (|201in . the stellar mass of DGl is in good agree- 
ment with galaxies at similar halo masses, as observed by 
THINGS. If high-temperature metal- line cooling had been 
added with eSN held constant, nDGl would have overpro- 
duced stars for galaxies in a comparable halo mass range. 
However, the enhanced feedback seems to have steepened 
the spatial power spectrum of the cold ISM of nDGl rela- 
tive to DGl, making it less consistent with the power spec- 
trum observed for the SMC. It is unclear, however, how the 
power spectrum varies with the instantaneous SFR and if 
this result holds across time. 

Capturing all the relevant ISM physics necessary to 
recover the full spectrum of turbulence sources at pc and 
sub-pc scales remains an outstanding challenge. Despite 
these limitations, the simulated dwarf galaxies presented 
here have been shown to possess bulk characteristics con- 
sistent with those observed in nature, including adherence 
to scaling relations such as the size-luminosit y, size- velocity, 
and luminosity- velocity IjBrooks et al.ll2011^ . Additionally, 
the star formation and feedback prescription used in these 
simulations has been shown to result in a realistic mass- 
metallicity relationship as a function of time, and consume 
gas at a rate that reproduces the incidence rate and metal- 
licities of both QS O-Damped Lyman Alpha (DLA) and 
GRB-DLA svstems (iBrooks et al.ll2007l : iPontzen et al.ll2008l : 



That is, although the microphysics of the ISM cannot be 
fully captured at the force resolutions that must be used 
currently in cosmological simulations, this does not largely 
impact the bulk macrophysics such as the rotation curves 
(stellar and dark matter mass profiles), angular momentum 
content, etc. On the other hand, we have seen that higher 
resolutions and adoption of more realistic physics for star 
formation leads to simulated galaxies that bett er reproduce 
the pro perties of observed g al axies (e. g.. iBooth et al 



2007 
2008: 



Robertson fc Kravtsovl 20081: Tasker fc Brvan 



Saitoh et all I2OO8I : ICeverino fc KlvpinI |2009| : 



Governato et al.l I2OI0I ). The work presented here high 



IPontzen fc PettinilbOOa T 

Hence, it is clear that our simulations remain ex- 
tremely successful in recovering many of the global optical 
and dynamical properties of realistic bulgeless dwarfs. 



lights paths for future improvement in the implementation 
of ISM physics in cosmological simulations, and provides 
useful tests for reassessment once metal-dependent H2 
cooling and star formation has been added to Gasoline 
and other cosmological simulation codes. 
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